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Abstract 

It is argued that large scale cosmic magnetic field could be generated in the 
primeval plasma slightly before hydrogen recombination. Non-zero vorticity, 
necessary for that, might be created by the photon diffusion in the second 
order in the temperature fluctuations. The resulting seed fields at galactic 
scale would be only 4 orders of magnitude smaller than the observed ones and 
with a mild galactic dynamo amplifying the seed fields by the factor ~ 10 4 an 
existence of coherent magnetic fields in galaxies may be explained. 

PA CS: 95. 30. Qd, 98. 62. En 
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1 Introduction 

Existence of magnetic fields on astronomically large scales remains one of unsolved 
cosmological mysteries. It is known from observations that there are magnetic fields 
in galaxies with the strength about micro-gauss which are homogeneous over galactic 
size I gai ~ ( a few) kpc, and maybe, even more puzzling, intergalactic magnetic fields 
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with the strength approximately three orders of magnitude weaker but homogeneous 
on much larger intergalactic scale of hundreds kpc (for review see e.g. refs. [1]). 
Though in stellar processes much stronger magnetic fields can be generated, their 
coherence scale is negligible in comparison with l ga i. There are two competing ideas 
of explanation of the phenomenon [1, 2]. The first one is based on traditional astro- 
physics and pursues the possibility that strong magnetic fields generated in multiple 
stellar catastrophes could be ejected into interstellar space and by field line recon- 
nection might create magnetic fields coherent on galactic size (for a recent review see 
ref. [3]). However there are no compelling quantitative arguments in favor of this 
hypothesis. Moreover, the energy density of galactic magnetic fields are of the same 
magnitude as the energy density of cosmic microwave background radiation (CMBR) 
so its magnitude is about 1CT 10 of the total galactic mass/energy. Such a huge con- 
tribution is difficult to explain by conventional mechanisms. It is even more difficult 
to explain in this way intergalactic magnetic fields, if they exist. 

Other possible mechanisms of generation of galactic magnetic fields are based on 
physical processes in the early universe - for the review see e.g. refs. [2, 4]. One can 
also find an extensive list of literature in ref. [5]. Basically there are three different 
mechanisms discussed: 

1. Breaking of conformal invariance of electromagnetic interaction at inflationary 
stage. The latter could be realized either through new non-minimal (and possibly 
non gauge invariant) coupling of electromagnetic field to curvature [6], or in dilaton 
electrodynamics [7], or by conformal anomaly in the trace of the stress tensor induced 
by quantum corrections to Maxwell electrodynamics [8]. 

2. First order phase transitions in the early universe [9] producing bubbles of new 
phase inside the old one. A different mechanism but also related to phase transitions 
is connected with topological defects, in particular, cosmic strings [10]. A recent 
discussion and a model of generation of large scale magnetic field can be found in 
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ref. [11]. 

3. Creation of stochastic inhomogeneities in cosmological charge asymmetry, either 
electric [12], or e.g. leptonic [13] at large scales which produce turbulent electric 
currents and, in turn, magnetic fields. 

In this work we will consider a new mechanism somewhat similar to those men- 
tioned in point 3 above but in contrast to them this mechanism does not demand any 
new physics and could be realized at relatively late stage of the universe evolution, 
namely, at red-shifts z ~ 10 3 — 10 4 , near hydrogen recombination. The mechanism 
suggested here is intermediate in time between the early universe ones and the as- 
trophysical mechanism which took place practically in contemporary universe after 
galaxies were formed and stellar explosions took place. 

The basic features of the suggested model of magnetic field generation are the 
following. We show that, despite low Reynolds number (3), non-zero, though small, 
vorticity can be generated in the cosmic electron-photon fluid due to inhomogeneity 
of the latter and especially due to different spectrum of inhomogeneities of electronic 
and photonic components. Such difference could be created at relatively late stage of 
cosmological evolution (near hydrogen recombination) even from initially adiabatic 
density perturbations. The motion of cosmic plasma would create electric currents 
because of different velocities of electrons and protons and since the vorticity of the 
motion would be non-vanishing the same would be true for the currents, Vx J^O. 
Such currents are known to generate magnetic fields. An attractive feature of the pro- 
posed mechanism is that no new physics has to be invoked. The mechanism operates 
in the standard cosmological model in the frameworks of the usual Maxwell electro- 
dynamics. In what follows we will estimate the magnitude of the generated magnetic 
field and show that it can be strong enough and have sufficiently large coherence scale 
so that after dynamo amplification (for the review of the latter see refs. [14]) it can 
explain the observed fields in the galaxies and, possibly, even intergalactic ones. 
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2 Hydrodynamics of cosmic plasma at recombina- 
tion 



Fluid motion under pressure forces is governed by the hydrodynamical equation (see 
e.g. the book [15]): 
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p (d t Vi + v k d k Vi) = -d^ + d k 



V (d k Vi + diV k - -SikdjV^ + di ((djVj) 



(1) 



where v is the velocity of the fluid element, p and p respectively are the energy 
and pressure densities of the fluid, and rj and ( are the first and second viscosity 
coefficients. In the case of constant viscosity coefficients this equation is reduced to 
the well known Navier-Stokes equation. The coefficient rj is related to the mean free 
path of particles in fluid as 

rj/p = v = l f (2) 

In what follows we disregard second viscosity (. 

The character of the solution to eq. (1) crucially depends upon the value of the 
Reynolds number 

ft = £ (3) 

where A is the wavelength of the velocity perturbation. If R ^> 1, then the fluid motion 
would become turbulent and non-zero vorticity would be created by spontaneously 
generated turbulent eddies. In the opposite case of low R the motion is smooth and in 
the case of incompressible fluid with homogeneous p and viscosity coefficients r\ and £ 
the fluid velocity would have potential character with vanishing vorticity, V x v = 0. 
However, the above assumptions of homogeneity is not precise and some vorticity can 
be generated even with low R. 

Let us consider now the cosmological epoch before the hydrogen recombination, 
when the plasma consists of three components: photons, electrons, and baryons 



mainly constituted by protons (for simplicity, we shall neglect small amount of 4 He 
nuclei). For temperatures T > 1 eV the energy density of (e,p, 7)-plasma is domi- 
nated by photons. Indeed, at the present day the CMBR energy density is p 7 ~ 0.26 
eV/cm 3 , i.e. approximately 10 3 smaller than the energy density of baryonic matter pb, 
and thus p b and p 7 become comparable at red-shift z 10 3 or at T « 0.23 eV. 1 The 
contribution of dark matter (DM) into energy density is not important for plasma 
hydrodynamics as far as the DM does not interact with plasma. For the moment, we 
also neglect cosmological expansion considering characteristic times smaller than the 
Hubble time H^ 1 at the appropriate epoch: 

= 3/2 27 kpC - , (4) 

where T e y = (T/l eV). At the present scale, it corresponds to (1 + z)H~ x = 
110T~^^ 2 [T e y + 0.76] -1 / 2 Mpc. This result is obtained for cosmological relativistic 
matter consisting from photons and neutrinos (contributing 68% with respect to pho- 
tons); the cosmological constant is not essential. For labeling the cosmological epoch 
we will interchangeably use the temperature T or red-shift 1 + z = 4260 T e y. 

As far as the cosmological plasma was dominated by photons, the viscosity coef- 
ficient v is determined by the photon mean free path: 

where <jt = 871a 2 /3m 2 = 6.65 x 10 -25 cm 2 is the Thomson cross-section, n e = /3n 7 = 
6 • 10~ 10 • 0.24T 3 is the electron number density, and X e is a fraction of the free 
electrons: X e (z) is practically 1 for z > 1500, and sharply decreases for smaller z's, 

reaching values ~ 10 -5 at z < 1000. For T > 1 eV the mean free path (5) is much 

1 Throughout this paper, we use the following values of the cosmological parameters: the CMBR 
temperature T 7 °' = 2.725 K, the baryon-to-photon ratio (3 — ns/n^ = 6x 1CP 10 , the Hubble constant 
h = 0.71, the matter density of the universe Q m h 2 = 0.135, and the baryon density Qbh 2 — 0.0224. 
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smaller than the Hubble horizon: 



Z 7 1.1 x 10 



-3 



1 + 0.76T r 



_i-|i/2 



cV 



(6) 



H- 1 X e (T)T cV 

This ratio becomes comparable to unity only at temperatures T ~ 0.3 eV, when 
X e C 1. At the present day scale the photon mean free path (5) would be lj® 
130 kpcX-^v 2 - 

Let us first estimate the Reynolds number of the fluid motion created by the pres- 
sure gradient in eq. ( 1 ) . To this end we assume that the liquid is quasi incompressible 
and homogeneous, so that the second term in the r.h.s. of this equation can be ne- 
glected. This is approximately correct and the obtained magnitude of fluid velocity 
is sufficiently accurate. In this approximation eq. (1) reduces to a much simpler one: 

9 t v+(vV)v-(/Av = -- (7) 

P 

A comment worth making at this stage. The complete system of equations in- 
cludes also continuity equation which connects the time variation of energy density 
with the hydrodynamical flux (see below eq. (19) and the Poisson equation for grav- 
itational potential induced by density inhomogeneities. We will however neglect the 
gravitational force and the back reaction of the fluid motion on the density pertur- 
bation. This approximation would give a reasonable estimate of the fluid velocity for 
the time intervals when acoustic oscillations are not yet developed, i.e. for t < X/v s , 
where A is the wave length of the perturbation and v s is the speed of sound (in the 
case under consideration v\ — 1/3). In fact the wave length should be larger than 
the photon mean free path, to avoid diffusion damping, and the characteristic time 
interval, as we see below, should be somewhat larger than A. So we may hope that 
our estimates of the velocity are reasonable enough. Neglecting gravitational forces, 
especially those induced by dark matter would result in a smaller magnitude of the 
fluid velocity, so the real effect should be somewhat larger. 
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For small velocities (or sufficiently small wavelengths) we may neglect the second 
term in the l.h.s. with respect to the third one. In this approximation the equa- 
tion becomes linear and can be easily solved for the Fourier transformed quantities. 
Assuming that the parameters are time-independent (though it is not necessary) we 
obtain: 

ik r i 
v k = $k [l - exp(-uk 2 t)\ (8) 

where 5^ = {Sp/p)k is the Fourier transform of relative density perturbations, Sp/p; 
its natural value is ~ 10~ 4 , though it might be much larger at small scales. The 
coefficient 1/3 comes from equation of state of relativistic fluid, p = p/3. 
Therefore, for the Reynolds number we obtain: 

5k 



3 (k v) 2 



1 - exp (-z/£; 2 t)] . (9) 



If Sk is weakly dependent on k, then R k is a monotonically rising function of the 
wavelength A = 2n/k. For t <C A 2 /V it takes the value 

R k = ^$k (10) 

Therefore, considering the perturbation with the wavelength A = 1-njk which 
enters the horizon at the time t ~ H^ 1 , we find that the maximum value of the 
Reynolds number is given by i?™ ax (if -1 /3Z 7 )5fc. For later times, the Reynolds 
number for the comoving wavelength decreases approximately as oc T until the density 
perturbation is completely damped by the photon diffusion. Thus, for T ^> 1 eV, 
using eq. (6), we obtain i?™ ax ~ 300<5fcT c v, and so for the development of turbulence 
one needs SkT e y > 0.1. This condition can be satisfied either at high temperatures, 
T > 100 eV, or for large density perturbations, 5 k 3> 10~ 4 . This is not the case 
for the adiabatic fluctuations with a nearly flat spectrum, as is generically predicted 
in inflationary scenarios. However, one cannot exclude a situation that at smaller 



wavelengths, corresponding to the present day scales < 100 kpc, there are substantial 
isocurvature fluctuations, e.g. the baryon density perturbations (5p/p)b are larger 
than the dark matter fluctuations (Sp/p)j)M- 

The baryon density perturbation of the wavelength A, after many oscillations, will 
be damped by the Silk effect [16] at the effective timescale t c s ~ A 2 / v. However, be- 
cause of a large Reynolds number, it could happen that during this time the photon 
diffusion dragging electrons would also produce non-zero vorticity. 2 In this case the 
generation of magnetic fields might be strongly amplified. We shall keep in mind an 
interesting possibility of large 5k at smaller scales, however now we shall mostly con- 
centrate on a more plausible possibility of small Reynolds numbers and non-turbulent 
flow of cosmic fluid in T ~ 1 eV range of temperatures. 

At first sight vorticity in laminar fluid motion is not generated in the approxima- 
tion given by eq. (7). However, this is not quite so. If both p and p are different 
functions of space points, pressure gradient may create motion with non-vanishing 
V x v. Indeed, from eq. (7) follows 



where f] = V x v and we assume that velocity is small so that the term quadratic 
in v was neglected. If the r.h.s. is non-vanishing, then ft would be non-zero too. 
However usually pressure density is proportional to the energy density, p = wp, with 
a constant coefficient w and hence V x (Vp/p) = 0. 

In the case under consideration, electrons are strongly coupled to photons. Their 
effective mean free path is 



2 The case of large Reynolds numbers, that might be created at earlier epoch, T ~ 1 MeV, due 
to large leptonic asymmetries, was considered in rcf. [13]. 




(11) 




(12) 



a Th n 7 



8 



The factor in the numerator takes into account large thermal momentum of electrons, 
(p) ~ y/Sm^T, so that they need many collisions, with momentum transfer Ap = T, 
for a significant change of their momentum. Still despite the presence of this factor, 
the mean free path of electrons is much smaller than that of photons simply because 
there are more than billion photons per a single electron in the plasma. So electrons 
are practically frozen in the plasma, while photons may diffuse to much larger dis- 
tance. One should keep in mind that the plasma must be locally electrically neutral 
and thus redistribution of electron inhomogeneities is accompanied by the same re- 
distribution of baryonic ones. This makes the mean free path of charged species of 
matter even smaller. On the other hand, a motion of homogeneous component of 
electron number density does not require dragging of baryons. 

For the wavelengths larger than mean free path of charged particles we may con- 
sider plasma as a fluid where charged particles are strongly coupled to photons. So 
for the estimate of the rotational velocity we will use the hydrodynamic equation 
(11). As we have already mentioned, vorticity, V x v could be non-zero only if the 
fluid is non-homogeneous and strictly speaking we should use eq. (1) with all param- 
eters depending upon space points. However, we believe that for a simple estimate 
eq. (11) may be sufficient. Since before recombination the interaction rate between 
radiation and charged particles is very high, plasma/liquid should be in local thermal 
equilibrium with all constituents having the same temperature T(x). If T would be 
the only parameter which determines the state of the medium, then vorticity would 
not be generated because we would have in our disposal only VT and it is impossible 
to construct non-vanishing V x v from the gradient of only one scalar function. How- 
ever, distributions of charged particles depend upon one more function, their chemical 
potential: 

(13) 



/ = exp 



E 



T(x) 
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where the dimensionless chemical potential £ can be readily expressed through particle 
number density n e pa ub = f3(x)n~, with (5{x) = 6 • 10~ 10 + 8f3{x): 



Hence we will find that the source term in the vorticity equation (11) is equal to: 



An essential feature here is that the spatial distribution of charged particles does 
not repeat the distribution of photons and hence the vectors Vp 7 and V/3 are not 
collinear. This could occur if for lambda corresponding to subgalactic scales, there 
exist baryon isocurvature fluctuations and thus p(x) and (3(x) have different profiles. 
As we have mentioned above, different mean free paths of photons and charged par- 
ticles would maintain such non-collinearity of the order of unity at the scales A ~ Z 7 . 
Moreover, even in the case of adiabatic perturbations a shift in the distribution of 
photons and charged particles could also be created because of acoustic oscillations 
that proceeded with different phases of radiation and matter densities. At the scales 
A < Z 7 perturbations in the the plasma temperature would be erased by the diffusion 
damping [16], while for A ^> Z 7 diffusion processes are not efficient and one would 
expect self-similar perturbation leading to collinearity of Vp 7 and V/3. On the other 
hand, when A entered under horizon acoustic oscillations begun which destroyed the 
self-similarity. Thus the expected wavelengths of vorticity perturbations should be 
between Z 7 < A < if -1 . 

We assume that there is no additional suppression of the source term (15) and 
by the order of magnitude its amplitude corresponding to the wavelength A can be 
evaluated as 



£(x) = In P(x) + const 



(14) 




(15) 




(16) 
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Taking 5T/T ~ 5(3/(3 ~ 3 • 1(T 5 we obtain S ~ lO" 8 ^ 1 . 

Now equation (11) can be solved in the same way as eq. (7) and we find: 



Vorticity could also be generated even if perturbations in plasma are determined 
by a single scalar function, for example, by T(t, x) because if might be proportional 
to the product <9jT(t,x) djT(t' ', x). These two gradients generally are not collinear if 
taken at different time moments t and t'. To see that, let us start from the Boltzmann 
equation for the distribution function f(t, x, E, p) of photons: 



where V = p/E is the particle velocity (not to be confused with the velocity v of 
macroscopic motion of the medium), for photons V = 1, while ti <C 1, £ and p are 
respectively the particle energy and spatial momentum, H is the universe expansion 
rate, F is an external force acting on particles in question (the latter is assumed to be 
absent), and I [f a , ...] is the collision integral depending on the distributions f a of 
all participating particles (for the definition of the I co n see e.g. eq. (43) of review [17]). 

At temperatures in eV-range only the Thomson scattering of photons on electrons 
is essential, so the collision integral is dominated by the elastic term. Integrating both 
parts of eq. (18) over d 3 p/ (2n) 3 we arrive to the continuity equation: 




(17) 



with Z 7 given by eq. (5). 




(18) 



n(x) + VJ = 



(19) 



where J is the photon flux given by 




(20) 
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and v is the average macroscopic velocity of the photon plasma. Using the standard 
arguments one can derive from eq. (19) the diffusion equation: 

n = DAn (21) 

where D Z 7 /3 is the diffusion coefficient. We will use this equation below to 
determine time evolution of the photon temperature T. 

If the elastic reaction rate T e i = a^n e X e = l// 7 is sufficiently large, local thermal 
equilibrium would be established and the photon distribution would be approximately 
given by 

/«/ = exp(- J E?/T + (22) 

where the temperature and effective chemical potential could be functions of time 
and space coordinates: T = T(£, x) and £ = £(£, x), and the photon mean free path 
is given by eqs. (5,6). Evidently fo annihilates the collision integral. We can find 
correction to this distribution, f — fo + fi, substituting this expression into kinetic 
equation (18) and approximating the collision integral in the usual way as — r e i/i: 

(K + T cl ) h = -Kf (23) 

where K is the differential operator, K = d t + (W). The solution of this equation 
is straightforward: 

r-t 



fi (t, x, E, V) = - / dn exp - / dr 2 T el (t - r 2 , x - Vr 2 ) 
jo L Jo 



■fl7o(f-Ti,x-VTi)(24) 



Using this result we can calculate the average macroscopic velocity of the plasma. The 
calculations are especially simple if elastic scattering rate is high and the integrals 
are dominated by small values of t±. In this case we obtain: 

JrfW l( t,x,E,V) 
V ' {t ' } ~ Sd?pf (t,x,E) (25) 
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and its vorticity, Q{ = eijkdjVk'- 



( 




T 



) 



( 



did t T 



T 



) 



(26) 



To estimate time derivatives of the temperature we will use the diffusion equa- 
tion (21), from which we find d t T = DAT and finally obtain for vorticity at the scale 



Since the photon diffusion erases temperature fluctuations at the scales A < Z 7 the 
vorticity has the maximum near A ~ l T This magnitude of vorticity is considerably 
larger than found previously (17) and we will rely on it in the estimates of magnetic 
field presented in the following section. 

To avoid confusion let us mention that at the moment when perturbation with a 
given wave length enters the horizon, A 3> Z 7 . Later (A// 7 ) scales as ~ a(t)~ 2 and till A 
remains larger than Z 7 the amplitude of perturbations does not decrease significantly 
and only when A < / 7 the photon diffusion damps density perturbations. 

3 Electrodynamics of cosmic plasma with helical 
flows and magnetic field generation 

Generation and evolution of magnetic field strongly depends upon the electric conduc- 
tivity of plasma which can be estimated as follows. Equation of motion of a charged 
particle in external electric field E has the usual form 



So the drift velocity gained during the time At is equal to Vd = eEAt/m. The charged 
particles (electrons) keep on to be accelerated approximately during time between 



A: 




(27) 



mV d = eE 



(28) 
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collisions, At = l e /V T where V T = ^JST/m is the thermal velocity and the electron 
mean free path l e is given by eq. (12). After collision the particle loses (forgets) its 
previous velocity and the process repeats. This is true for sufficiently weak fields when 
the drift velocity is small in comparison with thermal velocity, otherwise run-away 
charge carriers would be produced and the conductivity would be much larger. 

The induced electric current is J = en e V, where n e is the number density of 
charge carriers. Comparing this with the definition of current conductivity k — J/E 
we find: 

2a n 7 T y J 

The conductivity is very high, so the generation of magnetic field by the source 
currents, created by the cosmological inhomogeneities, is governed by the well know 
equation of magnetic hydrodynamics: 

d t B = V x (v x B) + - V x J (30) 

K 

The electric current J induced by the relaxation of the density inhomogeneities would 
contain two components: electronic and protonic. However the first one is surely 
dominant because it is much easier to drift electrons than heavier protons. This is 
why a non-zero current can be induced. Of course motion of electrons should not 
produce any excess of electric charge but this can be realized because the current is 
created in the dominant homogeneous part of the charge particle distribution. 
The solution of eq. (30) can be roughly written as 

B ~ fdtJ^-) e 2 ™' 1 ^ (31) 



A K 

The exponential factor under the integral, which presents pregalactic dynamo effect 
is normally rather weak, the exponent is about 

^ 500T oV (f)^ (32) 
14 



where eqs. (8) and (6) have been used and the density fluctuations were taken as 
5 X = 4{5T/T) X . With {ST/T) X ~ 3 • KT 5 and T ~ eV the exponent is about 
0.015 T e y < 1 and pregalactic dynamo is not important. It may be significant if the 
density perturbations at relatively small scales are much larger than their accepted 
canonical values. Exponential enhancement can be large at higher temperatures, T > 
10 3 eV, which correspond to the present day scales below kpc. Equipartition between 
magnetic field and CMBR can be expected on such scales and after 'Brownian" line 
reconnection and relatively mild galactic dynamo the observed galactic magnetic fields 
may be created. 

An estimate of magnetic field without pregalactic dynamo enhancement can be 
easily done if the helical source current is known, V x J = en e fl. With ft given by 
eq. (27) and B by eq. (31) we obtain: 

= 0.24- 10 W (1) (t) 3 W T 3, (33, 

where we took the wavelength equal to the photon mean free path, A = Z 7 . 

If we take into account that linear compression of pregalactic medium in the 
process of galaxy formation is approximately r ~ 10 2 , the seed field in a galaxy 
after its formation would be t 2 Bq, i.e. 4 orders of magnitude larger than that given 
by eq. (33) and, for T = 1 eV, a relatively mild galactic dynamo, about 10 4 , is 
necessary to obtain the observed galactic magnetic field of a few micro-Gauss at the 
scale Ib ~ (100/r) kpc = 1 kpc. The seed magnetic fields formed earlier (at higher T) 
would have larger magnitude (~ T 3 ) but their characteristic scale would be smaller by 
factor 1/T 2 . Chaotic line reconnection could create magnetic field at larger, galactic 
scale I ga i, but the magnitude of this field would be suppressed by Brownian motion 
law - it would drop by the factor (iB/lgai) 3 ^ 2 - It is interesting that according to our 
results all scales give comparable contributions at l ga \. This effect may lead to an 
enhancement of the field but it is difficult to evaluate the latter. Let us also note that 
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magnetic fields generated by the discussed mechanism at the cluster scale, 10 Mpc, 
should be not larger than 1CT 8 /iG if no additional amplification took place. 

4 Discussion and conclusion 

We have shown in this work that photon diffusion produced by temperature fluc- 
tuations of CMBR could create vortical flows in the primeval plasma in the second 
order, i.e. proportional to (ST/T) 2 . Because of different mobility of electrons and 
protons (or helium ions) in the plasma the helical macroscopic motion of the latter 
would create helical electric currents which would in turn generate magnetic fields. 
The amplitude of such fields generated at the moment when plasma temperature was 
equal to T is given by eq. (33). The characteristic scale of the field should be equal 
to the photon mean free path in the plasma at temperature T which at the present 
day is about 100 kpc (eV/T) 2 . This scale is quite close to the galactic size. Adiabatic 
compression in the process of galaxy formation by the factor r ps 100 would lead 
to an enhancement of magnetic field by r 2 = 10 4 . Thus to meet the observational 
data this seed field should be amplified by the galactic dynamo only by 4-5 orders 
of magnitude. Even if the magnitude of the seed field presented above is somewhat 
overestimated there is a huge reserve in the amplification by the galactic dynamo 
which could be as large as 15 orders of magnitude [14]. 

A nice feature of this mechanism is that it does not demand any new physics for the 
realization. Of course, density/temperature perturbations at large scales were, most 
probably, created during inflationary epoch and in this sense new physics is invoked 
but independently of the mechanism we know from observations that 5T/T ^ and, 
based on these data, we may estimate the seed magnetic fields using good old physics. 

It is interesting that, according to the existing indications to intergalactic magnetic 
fields, their strength is 3-4 orders of magnitude weaker than the strength of galactic 
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fields. This fact (if it is a fact) hints that intergalactic and galactic magnetic fields 
might have common origin, but galactic fields are larger due to mentioned above 
adiabatic compression which enhanced the field by 4 orders of magnitude. If it is 
so, the mechanism discussed in this work may be irrelevant because some galactic 
dynamo is needed to amplify the galactic seed field up to the observed magnitude. 
On the other hand, such dynamo seemingly does not operate on intergalactic scales. 

Larger density perturbations would be helpful for generation of larger magnetic 
field for which dynamo might be unnecessary. Though much bigger ST is not formally 
excluded at the scale about 100 kpc, but to have them at the level (ST/T) 2 ~ 10~ 4 
seems to be too much. A natural idea is to turn to a later stage, to onset of struc- 
ture formation when Sp/p becomes larger than 10" 2 . With such density perturbations 
strong enough magnetic fields may be generated without dynamo amplification. How- 
ever after recombination the number density of charge carriers drops roughly by 5 
orders of magnitude. Correspondingly Z 7 rises by the same amount and the strength of 
the seed field would be 5 orders of magnitude smaller if density perturbations and the 
temperature of formation remained the same. However both became very much dif- 
ferent. Density perturbations rose as scale-factor, (Sp/p) 2 ~ (T eq /T) 2 , where T eq ~ 1 
eV is the temperature when radiation domination changed into matter domination 
and density perturbations started to rise. Since, B/T 2 ~ T 3 , according to eq. (33), 
the net effect of going to smaller T is a decrease of B/T 2 which would be difficult to 
cure even by later reionization. Still, as argued in ref . [18] , magnetic field generation, 
driven by anisotropic and inhomogeneous radiation pressure (and in this sense similar 
to our mechanism) at the epoch of reionization, could end up with the field of about 
8 • 10~ 6 /iG. This result is 8 orders of magnitude larger than that found in the earlier 
papers [19] and quite close to ours (33), though these two mechanisms operated during 
very different cosmological epochs and were based on different physical phenomena. 

After this paper was sent to astro-ph we became aware of the work [20] " Gener- 
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ation of Cosmic Magnetic Fields at Recombination" where a much weaker effect was 
found. We think that the difference can be attributed to the following effects. We 
considered earlier period when the photon mean free path was much smaller than 
the horizon. It gives a factor about 10 3 in fluid velocity, eq.(8). Moreover, since in 
our case the electrons are tightly bound to photons the electron-photon fluid moves 
as a whole (while protons and ions are at rest) and the electric current induced by 
macroscopic motion/oscillations of plasma is noticeably larger. 

Acknowledgments 

We thank C. Hogan and P. Naselsky for critical comments. This work was partially 
supported by the MIUR research grant " Astroparticle Physics". 

References 

[I] P.P. Kronberg, Rep. Prog. Phys., 57 (1994) 325; 

R. Beck et al, Ann. Rev. Astron. Astrophys., 34 (1996) 155. 

[2] D. Grasso, H.R. Rubinstein, Phys. Repts., 348 (2001) 161. 

[3] P.L. Biermann, CF. Galea, to be published in Proc. Palermo Meeting, Sept. 
2002, Eds. N. G. Sanchez et al., astro-ph/0302168. 

[4] A.D. Dolgov, hep-ph/0110293; to be published in Collection of papers dedicated 
to 70th anniversary of S.G. Matinyan. 

[5] K. Dimopoulos, T. Prokopec, O. Tornkvist, A. C. Davis, Phys. Rev. D65 (2002) 
063505. 

[6] M.S. Turner, L.M. Widrow, Phys. Rev., D37 (1988) 2743. 



18 



[7] B. Ratra, Astrophys. J., 391 (1992) LI; 

D. Lemoine, M. Lemoine, Phys. Rev., D52 (1995) 1955; 

M. Gasperini, M. Giovannini, G. Veneziano, Phys. Rev. Lett., 75 (1995) 3796; 
M. Gasperini, Class. Quant. Grav., 17 (2000) Rl; 
M. Giovannini, hep-ph/0104214. 

[8] A.D. Dolgov, Phys. Rev., D48 (1993) 2499; 

A.D. Dolgov, ZhETF, 81 (1981) 417 [Sov. Phys. JETP, 54 (1981) 2223]; 

A.-C. Davis, K. Dimopoulos, T. Prokopec, O. Tornkvist Phys. Lett. B501 (2001) 

165; 

T. Prokopec, astro-ph/0106247. 

[9] C.H. Hogan, Phys. Rev. Lett., 51 (1983) 1488. 

[10] T. Vachaspati, A. Vilenkin, Phys. Rev. Lett., 67 (1991) 1057. 

[11] D. Boyanovsky, H.J. de Vega, M. Simionato, Phys.Rev. D67 (2003) 023502; hep- 
ph/0211022 (to be published in Phys.Rev. D). 

[12] A.D. Dolgov, J. Silk, Phys.Rev., D47 (1993) 3144. 

[13] A.D. Dolgov, D. Grasso, Phys. Rev. Lett. 88 (2002) 011301. 

[14] A.A. Ruzmaikin, A.M. Shukurov, D.D. Sokolov, "Magnetic Fields in Galaxies", 
Kluwer Academic Publishers, Dordrecht, 1988; 
R.M. Kulsrud, Ann. Rev. Astron. Astrophys. 37 (1999) 37; 
L. Malyshkin, R. Kulsrud, astro-ph/0202284. 

[15] L.D. Landau, E.M. Lifshits, "Gidrodinamika" (Hydrodynamics) (Moscow: 
Nauka, 1986) [Translated into English "Fluid Mechanics" (Oxford: Pergamon 
Press, 1987)]. 



19 



[16] J. Silk, Astrophys. J. 151 (1968) 459. 
[17] A.D. Dolgov, Phys. Rept. 370 (2002) 333. 

[18] M. Langer, J.-L. Puget, N. Aghanim, Phys. Rev. D67 (2003) 043505. 

[19] K. Subramanian, D. Narasimha, S.M. Chitre, Mon. Not. Roy. Astr. Soc. 271 
(1994) L15; 

N.Y. Gnedin, A. Ferrara, E.G. Zweibel, Astrophys. J. 539 (2000) 505. 
[20] C. Hogan, astro-ph/0005380. 



20 



